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ABSTRACT 


Specimens of the genus Orobdella Oka, 1895 from Korea, including various locations in the Korean Penin- 
sula, were identified as Orobdella tsushimensis Nakano, 2011. Phylogenetic analyses using mitochondrial cyto- 
chrome oxidase subunit 1 (COD, ND1, tRNA@’, tRNA™*, 12S rRNA, tRNA“, and 16S rRNA markers show 
that the newly collected specimens form a monophyletic group with the known O. tsushimensis specimens. 
The genetic distance of COI of these specimens was in the range 0.4—6.6%. These results confirm that the newly 
collected specimens belong to O. tsushimensis. This is the first record of the genus Orobdella from the 


Korean Peninsula. 
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Korean Peninsula 


INTRODUCTION 


The genus Orobdella Oka, 1895 is a terrestrial macrophag- 
ous leech taxon that is distributed in the Far East. At present 
11 species of this genus have been described: 10 from Japan 
(Oka, 1895; Richardson, 1975; Nakano, 2010, 2011a, 201 1b, 
2012a, 2012b, 2012c); and the other from Taiwan (Nakano 
and Lai, 2012). Additionally, Orobdella whitmani Oka, 1895 
has been reported from Primorsky Krai, Russia (Gilyarov et 
al., 1969), but the precise taxonomic status of the Russian 
specimens still remains uncertain (Nakano, 2012b). Yun 
(1997) provided a Korean leech list, but Orobdella leeches 
were not included in this list. Recently, quadrannulate Orob- 
della tsushimensis Nakano, 2011 has been recorded from 
Gageodo Island, which is located in the Yellow Sea, southwest 
of the Korean Peninsula (Nakano and Seo, 2012). However, 
Orobdella leeches have never been reported from the Korean 
Peninsula. 

Recently, additional specimens of the quadrannulate Orob- 
della were collected from various places in Korea, including 


locations in the Korean Peninsula. This is the first record of 
the genus Orobdella from the peninsula. The identification 
and a description of the newly collected specimens are herein 
presented. In addition, their phylogenetic relationships are 
reconstructed using mitochondrial cytochrome oxidase sub- 
unit 1 (COI), ND1, tRNA’, tRNA™", 12S rRNA, tRNA“, 
and 16S rRNA (tRNA@*-16S) sequence data. 


MATERIALS AND METHODS 


Leeches were collected from various places in Korea (Fig. 1). 
Altitudes and coordinates for localities were obtained using 
a Garmin eTrex GPS unit. 

For DNA extraction, botryoidal tissue was taken from the 
posterior part of the body around the caudal sucker of the 
specimens. All of the specimens were preserved in 70% etha- 
nol. Four measurements were taken: body length (BL) from 
the anterior margin of the oral sucker to the posterior mar- 
gin of the caudal sucker; maximum body width (BW); caudal 
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Fig. 1. Map showing the collection localities of the specimens examined in this study (nos. 1-11), Tsushima Island, Japan, and the 


type locality of Orobdella tsushimensis Nakano, 2011. 


sucker length (CL) from the anterior to the posterior margin 
of the sucker; caudal sucker width (CW) from the right mar- 
gin to the left margin of the sucker. Examination, dissection, 
and drawings of the specimens were accomplished using a 
stereoscopic microscope equipped with a drawing tube (Leica 
M125; Leica, Wetzlar, Germany). Specimens used in this 
study have been deposited in the Zoological Collection of 
Kyoto University (KUZ). 

The numbering convention is based on Moore (1927): body 
somites are denoted by Roman numerals and the annuli in 
each somite are given alphanumeric designations. 

The extraction of genomic DNA and DNA sequencing me- 
thods for COI and tRNA®-16S followed Nakano (2012b). 
For obtaining tRNA™ and ND1 (tRNA™'-ND1) sequence 
data, the procedure above was modified as follows: the primer 
sets for tRNA‘"-ND1 were LND300 and HND1932 (Light 
and Siddall, 1999) or HND1932-in (5'-AATGGAGCACGA 
TTAGTTTC-3’) (in this study); the PCR reaction mixture 
for tRNA-ND1 was heated to 94°C for 5 min, followed 
by 35 cycles of 94°C (10s), 45°C (20 s), and 72°C (39 s) and a 
final extension at 72°C for 6 min; the sequencing reaction 
mixture for tRNA™-ND1 was incubated at 96°C for 2 min, 
followed by 40 cycles of 96°C (10 s), 50°C (5 s), and 60°C 
(39 s). Obtained sequences have been deposited in GenBank. 
Sequences used for molecular phylogenetic analyses are 
listed in Table 1; 16 sequences of COI, 3 those of tRNA*- 
ND1 and 15 those of tRNA®-16S, which were published 
in the previous studies (Nakano, 2012a, 2012b; Nakano and 


88 


Lai, 2012; Nakano and Seo, 2012; Nakano and Gongalsky, 
2014), were obtained from GenBank. 

Sequences of mitochondrial COI and ND1 (tRNA™" region 
was removed from each sequence) were aligned by eye 
because there were no indels. Mitochondrial tRNA“*-16S 
sequences were aligned using MAFFT L-INS-I (Katoh et 
al., 2005) and then refined with GBLOCKS (Castresana, 
2000). The length of the aligned COI sequences was 1,266 
bp, that of ND1 was 579 bp, and that of tRNA@*-16S was 
818 bp (288 characters were eliminated). The concatenated 
sequences thus yielded a total of 2,663 bp positions. Pairwise 
comparisons of Kimura-2 parameter (K2P) distance (Kimura, 
1980) of the COI sequences were calculated using MEGAS 
(Tamura et al., 2011). 

Phylogenetic trees were constructed using maximum likeli- 
hood (ML) and Bayesian inference (BI) models. ML phylo- 
genies were calculated using TREEFINDER v October 
2008 (Jobb et al., 2004) with the tool package PHYLO- 
GEARS v 2.0 (Tanabe, 2008), and then nonparametric 
bootstrapping (Felsenstein, 1985) was conducted with 500 
replicates. The best-fit models for each partition were sele- 
cted based on the Akaike Information Criterion (Akaike, 
1974) by using KAKUSAN4 (Tanabe, 2011): for the Ist 
position of COI, TN93 with gamma distribution (+G) and 
proportion of invariant sites (+I); for the 2nd position of 
COI, TVM +I; for the 3rd position of COI, HK Y85+G; for 
the Ist position of ND1, GTR+G+I; for the 2nd position of 
ND1, J1+G; for the 3rd position of ND1, J2+G; and for 
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Fig. 2. Orobdella tsushimensis Nakano, 2011 from Seoul, KUZ Z616. A, Dorsal view; B, Ventral view. Scale bars: A, B=5 mm. 


tRNA“ -16S, GTR+G was selected. BI and Bayesian poste- 
rior probabilities (BPPs) were estimated using MRBAYES 
v 3.2 (Ronquist et al., 2012). The best-fit models for each 
partition were identified under the Bayesian information 
criterion (Schwarz, 1978) and also by using KAKUSAN4: 
for the 1st position of COI, GTR+G +I; for the 2nd position 
of COI, F81 +I; for the 3rd position of COI, HK Y85 +G; for 
the Ist position of ND1, GTR+G +I; for the 2nd position of 
ND1, HK Y85+G; for the 3rd position of ND1, HK Y85 +G; 
and for tRNA@*-16S, GTR +G was selected. Two indepen- 
dent runs of four Markov chains were conducted for 4 mil- 
lion generations, and the tree was sampled every 100 gene- 
rations. The parameter estimates and convergence were check- 
ed using TRACER v 1.5 (Rambaut and Drummond, 2009), 
and based on these results, the first 10,001 trees were dis- 
carded. 

Nodes with bootstrap (BS) values higher than 70% were 
considered as sufficiently resolved (Hillis and Bull, 1993). 
Nodes with BPPs higher than 95% were considered statisti- 
cally significant (Leaché and Reeder, 2002). 


SYSTEMATIC ACCOUNTS 


Order Arhynchobdellida Blanchard, 1894 
Suborder Erpobdelliformes Sawyer, 1986 
Family Gastrostomobdellidae Richardson, 1971 
**Genus Orobdella Oka, 1895 


** Orobdella tsushimensis Nakano, 2011 (Figs. 2, 3) 
Orobdella tsushimensis Nakano, 2011a: 41-46, figs. 2-4 

(near the Orobashi Bridge at Kechi Dam, alt. 69 m, 34° 

15'17’N, 129° 17'17"E, Tsushima, Tsushima Island, Naga- 


saki Prefecture, Japan, holotype, KUZ Z134); Nakano 
and Seo, 2012: 237-239, figs. 2-4. 


Material examined (see Fig. 1 for the locality number). 
Korea: 3 specimens, Seoul: Nowon-gu, Mt. Suraksan, alt. 
118 m, 37° 41'08"N, 127° 04'07"E (locality no. 1), 4 Jul 2012, 
Seo HY, KUZ Z615, Z616, dissected, and KUZ Z617; 2 
specimens, Jongno-gu, Mt. Bukhansan, alt. 447 m, 37°37’ 
28''N, 126°57'41"E (locality no. 2), 6 Jul 2012, Seo HY, 
KUZ Z621 and Z622, dissected; 2 specimens, Jeollanam- 
do: Sinan-gun, Gageodo Island, alt. 429 m, 34°04'49''N, 
125° 96'27"E, (locality no. 7), 10 Jul 2012, Seo HY, KUZ 
Z623 and Z624; 1 specimen, Wando-gun, Yeoseodo Island, 
alt. 46 m, 33°59'03"N, 126°55'19"E (locality no. 8), 20 Jul 
2012, Seo HY, KUZ Z628, dissected; 1 specimen, Jejudo 
Island: Jeju-si, alt. 384 m, 33°26’47"N, 126° 33'49"E (lo- 
cality no. 9), 17 Jun 2012, Seo HY, KUZ Z614; 1 specimen, 
Seogwipo-si, alt. 183 m, 33° 18'05"N, 126° 37'51"E (locality 
no. 11), 27 Sep 2012, Seo HY, KUZ Z629; 1 specimen, Jeju- 
si, Mt. Hanlasan, alt. 1,601 m, 33°21'56"N, 126° 30'48"E 
(locality no. 10), 2 Jul 2012, Kang SG, KUZ Z620, dissected; 
1 specimen, Gwangju Metropolitan City: Dong-gu, Mt. Mu- 
deungsan, alt. 533 m, 35°07'11''N, 126° 58'56"E (locality 
no. 6), 1 Jul 2012, Seo HY, KUZ Z618, dissected; 1 speci- 
men, Dong-gu, Mt. Mudeungsan, alt. 375 m, 35°07'21"N, 
126° 58'34"'E (locality no. 6), 24 Jun 2012, Seo HY, KUZ 
Z619; 1 specimen, Incheon Metropolitan City: Gyeyang-gu, 
Mt. Gyeyangsan, alt. 63 m, 37° 33'33"N, 126° 43'30"E (loca- 
lity no. 3), 11 Oct 2012, Seo HY, KUZ Z625, dissected; 1 
specimen, Gangwon-do: Wonju-si, Mt. Chiaksan, alt. 790 m, 
37° 18'07"'N, 128° 03'39"E (locality no. 4), 5 Oct 2012, Seo 
HY, KUZ Z626, dissected; 1 specimen, Daejeon Metropoli- 
tan City: Seo-gu, Mt. Changtaesan, alt. 245 m, 36° 11'38’'N, 
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Fig. 3. Orobdella tsushimensis Nakano, 2011 from Seoul, KUZ Z616. A, Dorsal view of somites I-VIII; B, Ventral view of somites |- 
VIII; C, Dorsal view of somites XXV-XXVII and caudal sucker; D, Ventral view of somites XXV-XXVII and caudal sucker; E, Ventral 
view of somites XI -XIII; F, Ventral view of gastropore and female gonopore; G, Ventral view of gastroporal duct; H, Dorsal view of 
epididymides and ventral nervous system; I, Dorsal view of male atrium showing position of ganglion XI. ac, atrial cornua; af, annu- 
lar furrow; an, anus; at, atrium; cp, crop; ed, ejaculatory duct; ep, epididymis; fg, female gonopore; gd, gastroporal duct; gp, gas- 
tropore; mg, male gonopore; np, nephridiopore; ph, pharynx; ts, testisac; vd, vas deferens. Scale bars: A-E, G-l=1 mm, F=0.5 
mm. 
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127° 20'19”E (locality no. 5), 5 Oct 2012, Seo HY, KUZ 
Z627, dissected. 

Description. Maximum BL 105.4 mm (KUZ Z622), BW 
11.9 mm (KUZ Z628), minimum BL 33.2 (KUZ Z626), BW 
4.3 mm (KUZ Z619). Caudal sucker ventral, ellipsoid, maxi- 
mum CL 3.9mm, CW 6.6 mm (KUZ Z622), minimum CL 
1.6mm (KUZ Z624), CW 2.0 mm (KUZ Z626). 

Somites VIII-XXV quadrannulate, somite XXVI triannu- 
late (Fig. 3A—E). Anus behind somite XX VII with no post- 
anal annulus (Fig. 3C). Clitellum in XI b5 to XIII a2. Eyes 
often in 3 pairs, first pair dorsally in posterior margin of II, 
second and third pairs dorsolaterally in posterior margin of 
V (al +a2) (Fig. 3A). Nephridiopores in 17 pairs in VII- 
XXIV (Fig. 3B, E). 

Pharynx reaching to XIII b5/b6-XIV b6 (Fig. 3G); each 
pharyngeal ridges reaching to crop straightly. Crop tubular, 
reaching to XIX b5/b6-XX al/a2. Gastropore conspicuous in 
XIII al/a2 (Fig. 3E, F). Gastroporal duct developed, bulbous, 
or bottle-shaped, slightly winding at junction with gastropore; 
joining with crop in XIII b6-XIV bó (Fig. 3G). 

Male gonopore in XI b6 (slightly anterior to middle of 
annulus), female gonopore inconspicuous in XIII al/a2, be- 
hind gastropore, gonopores separated by 1/2+5 annuli (Fig. 
3E, F). Testisacs multiple, in XVIII b5-XX b6 to XXIV al- 
XXVb6 (Fig. 3H). Paired epididymides in XV b5-XVIII al/a2 
to XVII b5-XX al/a2 (Fig. 3H). Ejaculatory ducts in XI b5 
to XV b5-XVIII al/a2; nearly straight in position anterior to 
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ovisacs (Fig. 3I); coiled in position posterior to ovisacs (Fig. 
3H). Pair of muscular atrial cornua, ovate in XI b5 and b6 
(Fig. 3D or fusiform in XI b6. Atrium short, muscular, glo- 
bular, in XI b5 and b6. One pair of ovisacs, thin-walled, glo- 
bular, in XIII a2 and b5. Oviducts thin-walled, both oviducts 
converging into common oviduct in XIII a2. 

Distribution. Known from mountainous regions of the Ko- 
rean Peninsula and its peripheral islands including Tsushi- 
ma Island, Japan (Nakano, 201 1a; Nakano and Seo, 2012). 
Phylogenetic relationships (see Fig. 1 for the locality num- 
ber). The ML tree with In L=—17134.56 (Fig. 4) was identi- 
cal in topology to the obtained BI tree (not shown). Mono- 
phyly of the genus Orobdella was confirmed (BS=100%, 
BPP=100%). Although O. tsushimensis formed a monophy- 
letic lineage with the clade that consists of ([O. esulcata Na- 
kano, 2010+0. mononoke Nakano, 2012]+[O. ketagalan Na- 
kano and Lai, 2012+ {0O. dolichopharynx Nakano, 2011 +0. 
shimadae Nakano, 2011}]), the phylogenetic position of O. 
tsushimensis was not statistically resolved (BS=45%, BPP= 
86%). 

Monophyly of O. tsushimensis was confirmed (BS=100%, 
BPP=100%). Within O. tsushimensis, three lineages were 
recognized and defined. Lineage A contained only KUZ Z626 
from Mt. Chiaksan, Gwangwon (locality no. 4). Lineage B 
consisted only of KUZ Z627 from Mt. Changtaesan, Dae- 
jeon (no. 5). Lineage C consisted of the other nine speci- 
mens. Lineages B and C formed a monophyletic lineage 
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Fig. 4. The maximum likelihood (ML) tree of 2,663 bp of cytochrome oxidase subunit 1 (COI), ND1, tRNA®S, tRNA™*, 12S rRNA, 
tRNA“, and 16S rRNA. The numbers associated with the nodes represent bootstrap values for ML(BS)/and Bayesian posterior 
probabilities (BPPs). BSs higher than 40% and/or BPPs higher than 95% are indicated. For the locality number, see Fig. 1 and Table 1. 
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(BS=91%, BPP=99%). Within Lineage C, a monophyletic 
lineage (BS=100%, BPP=100%) of three specimens, KUZ 
Z616 (no. 1), Z622 (no. 2) and Z625 (no. 3) from around 
Seoul, KUZ Z620 from Jejudo Island (no. 10) and a mono- 
phyletic lineage (BS=100%, BPP=100%) of two specimens, 
KUZ Z133 (no. 12) and Z134 (no. 13), from Tsushima Island 
formed a monophyletic lineage (BS=82%, BPP=99%). 
KUZ Z215 (no. 7), Z618 (no. 6) and Z628 (no. 8) from the 
southwestern part of Korea, including Gageodo and Yeoseo- 
do Islands, also formed a monophyletic lineage, but with 
low support values (BS=37%, BPP=58%). 

Genetic divergence. The COI K2P distance within O. tsu- 
shimensis was between 0.4-6.6% (mean=5.0%). 

Remarks. The present specimens newly collected from vari- 
ous places in Korea clearly belong to O. tsushimensis based 
on their possession of the following characteristics: somites 
VIII-XXV quadrannulate; gastropore conspicuous in XIII 
al/a2; gastroporal duct developed, bulbous; male gonopore 
in middle of XI b6, female gonopore in XIII al/a2, behind 
gastropore, gonopores separated by 1/2+5 annuli; paired epi- 
didymides in XV —XX, occupying more than 2 somites; atrial 
cornua developed, ovate or conical. One specimen, KUZ 
Z627, possesses tubular gastroporal duct and fusiform atrial 
cornua. However, these characteristics of this specimen seem- 
ed to be due to its immaturity, since this individual has unde- 
veloped sperm ducts and undetectable testisacs. In Nakano 
and Seo (2012), O. tsushimensis has been diagnosed as hav- 
ing epididymides in XVII-XIX (occupying approximately 2 
somites). However, the present specimens have their epidi- 
dymides in somites XV —XX (occupying 8-12 annuli, i.e., 2— 
3 somites). Therefore, the diagnosis in Nakano and Seo (2012) 
should be amended to take this difference into account. 

The molecular phylogenies obtained here showed that the 
specimens that were morphologically identified as O. tsushi- 
mensis form a well-supported monophyletic group. Addi- 
tionally, the COI sequence divergence within these speci- 
mens was between 0.4-6.6% (mean=5.0%). For compari- 
son, the intraspecific variation of COI sequences of O. koikei 
in Hokkaido, Japan was between 4.8-8.1% (mean=7.1%) 
(Nakano, 2012b). Thus, our molecular phylogenetic analy- 
ses and genetic distance data strongly support our morpho- 
logical identification result that the Orobdella leeches new- 
ly collected in this study belong to O. tsushimensis. 

The present molecular phylogenies showed that the speci- 
men from Mt. Chiaksan (Lineage A; locality no. 4) first 
diverged from the others. Then, the specimen from Mt. 
Changtaesan (Lineage B; no. 5) diverged from the other lee- 
ches. It is noteworthy that the specimen from Jejudo Island 
(no. 10) and two specimens from Tsushima Island (nos. 12, 
13) are monophyletic, which was statistically supported, with 
the three specimens from locations around Seoul (nos. 1-3), 
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not with the specimens from adjacent localities (i.e., nos. 6- 
8). This geographic pattern of the phylogenetic relationships 
suggests a complex biogeographical history of O. tsushi- 
mensis. However, the present phylogenies failed to reconstruct 
fully robust phylogenetic relationships within Lineage C. 
Additional specimens of O. tsushimensis as well as suitable 
genetic markers are still needed to clarify the precise phylo- 
genetic relationships and biogeographical history of this spe- 
cies. 
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